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The  function  of  a  data  acquisition  system  is  to  digitise  data  from  a  number  of  separate 
channels,  and  store  these  samples  for  later  analysis.  Sonar  systems  are  typically  multi¬ 
channel,  and  the  associated  data  acquisition  systems  are  required  to  handle  data  at  high 
data  rates.  The  high  data  rates  involved  in  digitising  bandpass  sonar  data  usually  mean 
that  the  only  viable  storage  methods  available  are  large  magnetic  tapes,  with  data  recorded 
on  high-speed  tape  recorders.  One  of  several  disadvantages  of  tms  type  of  storage  method 
is  that  the  types  have  to  be  replayed  through  special  decoders  to  obtain  the  data  in  a  form 
ready  for  analysis.  This  is  a  time  consuming  process,  delaying  the  analysis  and 
interpretation  of  the  data.  An  alternative  method  of  sampling  bandpass  sensor  data  which 
results  in  a  much  lower  data  rate  is  presented  in  this  report.  The  resulting  lower  data  rate 
allows  advantage  to  be  taken  of  the  recent  improvements  in  computer  storage  technology, 
enabling  a  data  acquisition  system  to  be  constructed  that  is  computer-based,  and  utilises 
computer  storage  devices  as  its  primary  storage  media.  The  advantage  of  this  storage 
method  is  that  the  data  is  in  computer  readable  form,  ready  for  immediate  analysis.  This 
work  was  carried  out  in  support  of  Navy  Task  NAV  92/401  and  an  application  of  the 
theory  of  the  AN/SQS  56  sonar  is  presented  in  the  Appendices. 
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1  Introduction 


The  function  of  a  data  acquisition  system  is  to  digitise  data  from  a  number  of  separate 
channels,  and  store  these  samples  for  later  analysis.  The  majority  of  active  sonar  systems 
have  a  large  number  of  channels,  resulting  in  high  overall  data  rates  when  the  signals 
from  these  channels  are  digitised  for  storage.  In  the  past,  the  only  storage  medium 
capable  of  handling  these  high  data  rates  was  large  magnetic  tape  used  with  high  speed 
tape  recorders  and  specialised  encoding  hardware.  This  approach  generates  recorded 
data  formats  which  require  further  processing  through  decoders  to  convert  them  to 
a  computer  readable  form,  making  it  costly  in  both  time  and  money.  However,  an 
alternative  approach  eliminating  some  of  the  drawbacks  of  this  method  is  now  possible, 
and  is  the  subject  of  this  report. 

Computer  storage  technology  has  reached  the  point  where  it  is  possible  to  store  in  excess 
of  one  gigabyte  of  data  on  a  single  disk,  with  sustained  data  rates  of  hundreds  of  kilobytes 
per  second.  With  this  level  of  technology,  it  is  now  feasible  to  construct  a  cost-effective 
data  acquisition  system  for  some  sonar  systems  where  the  digitised  data  is  stored  directly 
in  computer-readable  form.  Such  a  system  eliminates  the  need  for  pre-processing  of  the 
recorded  data  to  obtain  the  data  in  a  computer-readable  form. 

This  report  outlines  a  method  for  efficiently  digitising  bandpass  sonar  signals.  A  data 
acquisition  system  using  this  method  would  be  able  to  use  current  computer  storage 
technology  to  store  the  digitised  data  in  computer-readable  form.  Section  2  gives  a  brief 
introduction  to  bandpass  signals,  describing  how  they  can  be  represented  in  terms  of 
their  lowpass  complex  envelope,  and  how  bandpass  sampling  techniques  can  be  utilised 
to  sample  these  signals  at  a  rate  consistent  with  twice  the  signal  bandwidth.  Section  3 
outlines  a  method  of  obtaining  a  sampled  representation  of  the  complex  envelope  that 
is  only  applicable  to  a  certain  class  of  bandpass  signals. 

An  outline  of  how  these  techniques  can  be  applied  to  make  a  data  acquisition  system  for 
bandpass  sonar  data  which  is  efficient,  and  uses  computer-based  storage,  is  presented 
in  Section  4.  A  summary  is  given  in  Section  5.  An  application  of  the  data  acquisition 
method  detailed  in  this  report  to  the  AN/SQS-56  hull  mounted  sonar  is  described  in 
Appendices  C  through  D. 

It  is  important  to  note  that  the  mathematics  in  this  report  is  not  presented  rigorously, 
but  are  presented  in  such  a  way  as  to  enable  the  concepts  to  be  understood  easily.  A 
more  rigorous  mathematical  approach  to  the  techniques  outlined  can  be  found  in  the 
references  given  in  the  bibliography. 


2  Bandpass  Signals 

2.1  Complex  Envelope  of  Bandpass  Signals 

Consider  a  bandpass  signal  r(t)  of  bandwidth  W,  centred  at  frequency  f„  with  f0  »  W. 
Let  X(f)  denote  the  Fourier  transform  of  the  signal  x(t),  which  for  x(t)  real  is  defined 
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as  follows 
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X-(f)  and  X+(f)  respectively  denote  the  negative  and  positive  halves  of  the  spectrum 
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This  spectrum  is  shown  in  figure  1  below 
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Figure  1:  Spectrum  of  bandpass  signal 

It  is  well  known  [1],  [2]  that  a  bandpass  signal,  such  as  x(t)  above,  can  be  represented 
in  terms  of  its  in-phase  and  quadrature  components,  as  given  below 


x(t)  =  xi(t)  cos(27r  f0t )  —  x<j(f)  sin(2:r  fat)  (2) 

where  x/(t)  and  xg(t)  are  respectively  the  in-phase  and  quadrature  components  of  x{t). 
If  x(t)  has  bandwidth  W,  then  the  in-phase  and  quadrature  components  of  x(<)  will  each 
be  lowpass  signals  of  bandwidth  — . 

The  in-phase  signal  x/(t)  is  typically  obtained  by  modulating  the  signal  x(t)  by  a  carrier 
signal  of  frequency  f0  (2cos(27 rf0t)),  and  lowpass  filtering  to  remove  the  sum  frequency 
term.  Similarly  xq(<)  is  obtained  by  modulating  x(<)  with  a  carrier  90°  out  of  phase 
with  the  previous  carrier  (2sin(27r/0t)),  and  lowpass  filtering,  i.e. 

xi(t)  =  2cos(27r/0t)x(<)j^  (3) 

and 

xQ(t)  =  2sin(27r/0*)x(<)|  (4) 
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Figure  2:  Conventional  Quadrature  Demodulation 


where 


denotes  the  lowpass  filtering  operation.  The  conventional  method  of  obtaining 


the  in-phase  and  quadrature  signals  is  via  two  analog  multipliers  or  mixers  and  a  single 
90°  phase  shift,  as  shown  in  figure  2. 


The  in-phase  and  quadrature  components  of  the  bandpass  signal  can  be  combined  as 
follows,  to  form  a  complex-valued  signal  xc(t),  which  is  the  complex  envelope  of  the 
original  bandpass  signal. 


Xc(t)  =  Xi(t)  +jxQ(t)  (5) 

Hence  all  the  information  in  the  original  bandpass  signal  x(t )  is  contained  in  the  complex 
valued  lowpass  signal,  xc{t),  of  bandwidth 

2.2  Bandpass  Sampling  Techniques  for  Bandpass  Signals 

If  the  complex  signal  is  sampled  at  a  rate  consistent  with  its  bandwidth  of  y  (i.e 
sampled  at  W  samples  per  second)  then  the  complex  samples  provide  a  discrete-time 
representation  of  the  original  bandpass  signal,  with  a  sample  rate  consistent  with  the 
signal  bandwidth.  This  provides  a  much  more  efficient  sampling  method  (in  terms  of 
data  rate)  for  bandpass  signals  than  conventional  lowpass  sampling  techniques,  which 
require  a  sampling  rate  of  twice  the  highest  frequency  of  interest. 

The  lowpass  complex  signal  xc(t)  is  typically  sampled  by  forming  the  two  real  signals 
xi(t)  and  xg(f),  as  shown  in  figure  2,  and  then  sampling  each  of  these  signals  at  a  rate 
of  W  samples  per  second.  This  gives  a  combined  data  rate  of  2 W  real  samples/second 
or  W  complex  samples  per  second.  One  of  the  problems  of  sampling  xc{t)  using  this 
method  is  that  DC  biases  are  often  introduced  into  the  signals  x/(f)  and  XQ(t)  due  to 
the  non-ideal  nature  of  the  analog  mixing,  which  can  cause  problems  when  the  in-phase 
and  quadrature  signals  undergo  further  processing.  This  problem  is  often  overcome 
by  sampling  the  original  bandpass  signed  using  conventional  sampling  techniques,  and 
performing  the  quadrature  demodulation  digitally  as  shown  in  figure  3  below. 
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Figure  3:  Digital  complex  demodulation  of  a  bandpass  signal 

Other  bandpass  sampling  techniques  exist  for  efficient  sampling  of  bandpass  signals  such 
as  sampling  of  the  analytic  signal  (i.e.  the  sum  of  the  signed  plus  its  Hilbert  transform) 
and  second  order  sampling  (sampling  of  the  signal  and  its  quarter  wave  translation  at 
rates  consistent  with  the  signal  bandwidth).  More  details  on  these  techniques,  and  on 
quadrature  sampling  can  be  found  in  references  [1]  -  [6]. 


3  Quadrature  Sampling  -  An  Alternative  Method 

3.1  Integer  Band  Shifting 

For  bandpass  signals,  it  can  be  shown  that  if  a  suitable  relationship  exists  between  the 
signal  bandwidth  and  the  highest  signal  frequency,  then  a  sampling  frequency  can  be 
chosen  which  need  not  be  any  greater  than  twice  the  signal  bandwidth,  instead  of  twice 
the  highest  frequency  of  the  signal.  Consider  the  spectrum  of  the  bandpass  signal  shown 
in  figure  1.  Let  fy  denote  the  upper  band  edge  and  fy  denote  the  lower  band  edge. 


II 
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-f 

(6) 

,  ,  W 

h-fo~Y 

(7) 

< 

1 

II 

(8) 

where  f0  is  the  centre  frequency  of  the  band  and  W  the  bandwidth.  The  upper  and 
lower  band  edges,  fy  and  fy,  include  guard  bands  (i.e  filter  transition  bands)  so  they  in 
effect  denote  the  edges  of  the  stop  bands  of  the  bandpass  signal. 

Integer  band  shifting  is  achieved  by  intentionally  under-sampling  the  bandpass  signal 
with  respect  to  its  highest  frequency,  but  not  its  bandwidth.  Use  of  a  suitable  sampling 
frequency  results  in  the  positive  and  negative  halves  of  the  original  bandpass  spectrum. 
F+  and  respectively,  being  shifted  down  to  base  band  without  ambiguous  aliasing. 
There  are  two  possibilities  for  a  suitable  sampling  frequency.  The  first  is  where  the 
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sampling  frequency  f,  and  the  lower  and  upper  band  edges  ft  and  fc  are  related  as 
follows 


II  =  kf,  (»• 

fu=(^k  +  ^jfs,  where  k  =  1,2, 3,4, ...  (10) 


The  second  is  where  the  lower  and  upper  band  edges  /i  and  fc  relate  to  the  sampling 
frequency  /,  as  follows 


h  =  (*")/.  ,1U 

fc  =  kf,,  where  k  =  1,2,3. 4 _  (12) 

To  avoid  aliasing,  the  Nyquist  criterion  requires  that 

f,  >  2(fc  -  fL)  (13) 

Because  fc  and  fc  include  guard  bands,  equality  actually  applies  in  equation  ( 13),  giving 


f,  =  2 (fc  -  h) 


(14) 


iX(f) 
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Figure  4:  Integer  Band-shifting  of  Bandp^s  Signals 

Figure  4  shows  an  example  where  (11)  and  (12)  apply.  In  the  example  shown,  the  upper 
frequency  limit  fc  is  three  times  the  sampling  frequency  /,.  X (/)  and  S(f),  shown  in 
figures  4(a)  and  4(b)  are  the  spectra  of  the  original  bandpass  signal  and  the  sampling 


signal  respectively.  Figure  4(c)  is  the  spectrum  of  the  sampled  bandpass  signal,  which 
can  be  found  by  convolving  A .(f)  and  S{f).  Figure  4(c)  shows  that  the  spectrum  of  the 
sampled  signal  results  in  F.  and  F_  being  shifted  down  to  baseband  without  ambiguous 
aliasing. 

It  is  of  interest  :c  ^te  that  the  application  of  the  second  solution,  and  in  particular 
equation  12.  results  in  the  positive  half  of  the  original  spectrum  F+  being  shifted  to  the 
negative  half  of  the  resulting  spectrum,  and  F_  being  shifted  to  the  positive  half.  This 
doe?  not  have  any  undue  effect  on  the  further  processing  of  the  sampled  signal.  When 
(9)  and  (10)  are  used.  F+  and  F_  are  shifted  to  the  positive  and  negative  sides  of  the 
resulting  spectrum  respectively. 

Because  of  the  somewhat  restrictive  relationship  required  for  this  sampling  technique 
to  be  useable  (between  the  upper  and  lower  band  edges,  fi  and  fc .  and  the  sampling 
frequency.  fs.  as  shown  in  equations  (14)  and  either  (9)  and  (10).  or  (11)  and  (12) 
),  it  is  not  a  generally  applicable  technique  and  therefore  is  not  well  known  or  widely 
documented.  Some  details  of  this  technique  can  be  found  in  [7]  -  [9]. 


3.2  An  Alternative  Technique  for  Complex  Demodulation  us¬ 
ing  Integer  band  shifting 

If  the  nature  of  the  bandpass  signal  is  such  that  integer  band  shifting  on  sampling  is 
possible,  (as  described  in  section  3.1),  then  another  technique  [9]  exists  for  obtaining 
sampled  versions  of  the  in-phase  and  quadrature  signals.  This  alternative  technique 
avoids  the  analog  mixing  and  filtering  processes  (and  the  associated  problems)  required 
for  conventional  quadrature  demodulation  as  shown  in  figure  2.  and  also  avoids  the  high 
data  rates  associated  with  the  digital  technique  which  is  shown  in  figure  3. 

Consider  the  technique  outlined  in  figure  5. 
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Figure  5:  An  alternative  technique  for  complex  demodulation  using  integer  band  shifting 

The  spectrum  of  the  sampled  signal  x(nT)  in  this  figure  is  similar  to  that  shown  in 
figure  4(c).  The  spectrum  of  this  sampled  signal  can  be  defined  as  follows  ( using  notation 
similar  to  that  developed  in  Appendix  A) 

x(nT)  <=►  F-{j  — 2”  ±  2 mW^j  +  F+  +  —  ±  2mli'j  m  an  integer  (1-5) 
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where  <=>  denotes  the  Fourier  transform  operation.  F-  yf  —  5  *•  the  positive  half  of 

the  spectrum  of  the  sampled  signal,  and  F+  is  the  negative  half. 

The  in-phase  and  quadrature  bandpass  filters.  /i/(n)  and  hQ{n)  respectively,  are  both 
derived  from  the  same  lowpass  filter,  /i£.(n),  of  bandwidth  ~ .  via  the  following  frequency 
transformations. 


'W  / 

X  —  1  \  \ 

(16) 

hi(n)  =  hL(n)c os( 

,  2  V  ~  2  ) 

.  / 

'W  f 

V  -  1  \  \ 

(IT) 

hQ{n)  =  hi(n)  sinl 

,?r“  2  )) 

where  .V  is  the  order  of  the  lowpass  filter  hi(n). 

The  two  bandpass  filters  will  have  identical  amplitude  responses  but  will  be  90°  out  of 
phase.  To  simplify  the  analysis,  assume  that  the  lowpass  filter.  hi{n),  has  a  frequency 
response. 


HlU)  =  H(f  ±  m /,)  m  an  integer 
=  H(f±2m\V) 


(IS) 


since  f,  =  2H  .  and  where 


H(f)  =  { 


1. 

0  otherwise 


(19) 


The  frequency  response  of  the  two  bandpass  filters  will  be  as  follows 


Hi{f)  =  H^f ±2m\v'j  +  H^f  +  —  ±2mlVSJ,  m  an  integer  (20) 
Hq{/)  =  -]H{j  -  ±  2mW^j  +  jH  (f  +  —  ±  2mlV'^.  m  an  integer  (21) 

Consider  first  x(nT)  being  filtered  by  the  in-phase  bandpass  filter  k/(n).  The  resulting 
spectrum  (denoted  Xjf{/)  for  in-phase  filtered  version  of  x(nT),  xif{kT))  will  be 
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XiF{f)  =  ///(/)  *  (*-(/- y± 2 mw)  +  F+  (/  +  y  ±  2m W')  ) 

=  //(/-  y  ±  2miv)  F_  (/  -  ^  ±  2m^) 

/  W  \  (  W  \ 

+  H[f  +  —  ±2mW\F+[J  +  —  ±2mlVj 

The  signal  X[F(nT)  is  then  decimated  by  a  factor  2,  reducing  the  sample  rate  from  2 IT 
to  W  samples/sec.  This  results  in  the  portion  of  the  spectrum  that  covered  the  range 
— W  <  /  <  +W  being  replicated  at  intervals  of  IT,  instead  of  21V.  The  spectrum  of  the 
decimated  signal  xjf(nT')  (where  T'  =  2 T)  will  be 


(22) 


xIF(nV)  **#(/-  y  ±  mW^F.  (/  -  y  ±  mW') 

/  W  \  W  \ 

+  Hyf  — —  ±  mWjFM  ~  y  ±  mlVj  (23) 

Note  that  this  expression  shows  that 

H (/  -  y  ±  mw)F_  (/  -  y  ±  mw) 

and 

tf(/-y  ±mw)F+(/-y  ±mw) 

axe  aliased.  Finally  the  signal  xiF(nT')  is  multiplied  by  (—1)”,  which  in  the  frequency- 
domain  is  equivalent  to  shifting  the  spectrum  by  (or  r  radians).  Thus  the  spectrum 
of  the  signal  after  this  operation,  designated  xi(nT')  in  figure  5,  is 


•  • 


xi(nT')  Xi(f) 

Xi(f)  =  H(f  ±  mW)F-U  ±  mW)  +  H(f  ±  mW)F+{f  ±  mW)  (24) 


Ignoring  the  spectral  replication  at  intervals  of  ±mlT  which  exists  because  of  the  nature 
of  sampled  signals,  this  becomes 


X,(f)  =  mf)F.(f)  +  H(f)F+(f)  (25) 
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If  H(f )  is  an  ideal  lowpass  filter  of  bandwidth  y,  as  described  in  equation  (19),  then 
comparing  equation  (25)  with  equation  (A-7)  in  Appendix  A,  it  can  be  seen  that  xj(nT') 
is  a  sampled  version  of  x/(t),  the  in-phase  component  of  the  original  bandpass  signal. 


Following  a  similar  argument,  the  spectrum  of  xqF(tiT)  (designated  XQF(f)}  is 


XQF(f)  =  HqU)  X  (f.  (/  -  y  ±  2  mW^j  +  F+  (/  +  y  ±  2  mW^j 

=  -jH[f  -  -  ±  2mW)F-(f  -  y  ±  2m»'J 

+  ;//(/  +  y  ±  2mW'j  F+  (/  +  y  ±  2mwj  (26) 

xqF(tiT)  is  then  decimated  by  a  factor  of  2,  reducing  the  sample  rate  to  W  samples/sec. 

""Ms  gives  iQF(nr')  (where  T'  —  2T)  which  has  a  spectrum  as  follows 

xQF(nT')  -jH (/  -  y  ±  mwy _  (/  ~  y  ±  miv) 

+  jH{f  ~  f  ±  m^')F+  (/  "  Y  ±  mVr)  (2?) 

Finally  the  signal  XQF(nT')  is  multiplied  by  (-l)n  giving  the  signal  designated  as 
xq(tiT')  in  figure  5,  which  has  the  spectrum 

xQ(nT')  *  XQ(f) 

XqU)  =  ~]H{f  ±  mW)F.(f  ±  mW )  4-  jH(f  ±  mW)F+(f  ±  mil')  (28) 

Ignoring  the  spectral  replication  at  intervals  of  ±mW,  this  can  be  written  as  •  * 


•  • 


XQ(f)  =  -jHU)F-(f)  +  jH(f)F+(f)  (29) 

As  with  the  previous  case,  if  H(f)  is  eua  ideal  lowpass  filter  of  bandwidth  as  de¬ 
scribed  in  equation  (19),  then  comparing  equation  (29)  above  with  equation  (A-S)  in 
Appendix  A,  it  can  be  seen  that  xq(tiT')  is  a  sampled  version  of  XQ(t),  the  quadrature 
component  of  the  original  bandpass  signal. 

The  results  shown  in  equations  (25)  and  (29)  show  that  the  technique  outlined  in  figure  5 
is  a  valid  method  of  obtaining  sampled  versions  of  the  in-phase  and  quadrature  compo¬ 
nents  of  a  bandpass  signal,  provided  that  equations  (13)  and  either  (9)  and  (10),  or  (11) 
and  (12)  are  satisfied.  It  is  important  to  note  that  the  same  result  can  be  achieved  by 
removing  the  in-phase  bandpass  filter,  and  replacing  the  quadrature  bandpass  filter  with 
a  Hilbert  transformer  [4].  However  the  method  described  above  gives  superior  matching 
between  the  quadrature  channels. 
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4  An  Efficient  Data  Acquisition  System  for  Band¬ 
pass  Signals 


As  the  required  product  of  a  data  acquisition  system  is  computer- readable  data,  it  would 
be  highly  desirable  for  the  system  to  store  the  data  directly  in  a  computer-readable  form. 
In  the  past,  this  has  not  been  possible  due  to  the  high  data  rates  involved  and  the  limited 
data  rate  and  storage  capacities  of  computer  hardware.  However,  computer  technology 
has  reached  the  point  where  storage  capacity  in  excess  if  1  Gbyte  is  available  on  disk 
drives,  with  sustained  throughput  to  these  devices  in  the  order  of  hundreds  of  KBytes 
per  second.  If  a  sampling  scheme  is  used  where  the  bandpass  data  is  sampled  at  a  rate 
of  twice  its  bandwidth,  then  much  lower  data  rates  will  be  achieved.  Combining  such 
a  sampling  scheme  with  the  current  generation  of  computer  hardware  will  give  a  data 
acquisition  system  whose  stored  output  is  computer-readable,  and  ready  for  immediate 
analysis.  An  added  bonus  is  that  there  are  now  fewer  steps  in  the  process,  since  the 
decoding  step  is  redundant,  and  consequently  there  are  cost-savings  in  employing  this 
method. 

Consider  the  technique  of  complex  demodulation  outlined  in  section  3.2.  This  tech¬ 
nique  first  samples  the  bandpass  data  at  a  rate  of  twice  the  signal  bandwidth,  and 
then  processes  this  sampled  data  further  to  obtain  sampled  versions  of  the  in-phase  and 
quadrature  components  of  the  bandpass  signal.  If  this  technique  is  used  as  part  of  a  data 
acquisition  system,  then  the  data  could  be  stored  immediately  after  sampling,  with  the 
remainder  of  the  complex  demodulation  process  performed  at  some  later  stage.  As  this 
method  uses  a  sample  rate  of  twice  the  signal  bandwidth,  it  can  be  used  in  an  acquisition 
system  that  is  computer-based  and  uses  computer  storage  devices  to  store  the  results 
of  the  acquisition  process.  It  is  important  to  note  that  this  sampling  technique  is  only 
applicable  to  bandpass  signals  that  satisfy  equations  (13)  and  either  (9)  and  (10)  or  (11) 
and  (12)  (see  sections  3.1  and  3.2). 

A  possible  architecture  for  such  a  data  acquisition  system  is  shown  in  figure  6. 

In  this  system,  the  N  multi-channel  data  acquisition  cards  are  connected  to  a  computer 
bus  (e.g.  VME  bus).  Also  connected  to  this  bus  is  a  controlling  CPU.  which  has  multiple 
disk  drives  and  a  tape  drive  as  peripherals.  The  controlling  CPU  is  running  a  software 
process  that  controls  the  data  acquisition  cards.  As  the  data  is  captured  by  the  cards, 
the  CPU  reads  the  data  over  the  bus,  and  then  writes  the  data  to  a  file  on  one  of  the 
disk  drives.  This  process  continues  until  the  disk  drive  in  use  becomes  full,  at  which 
stage  the  file  writing  process  begins  writing  to  the  second  disk.  This  activity  continues 
until  all  the  disk  capacity  has  been  utilised.  At  this  point,  data  acquisition  will  cease, 
and  mother  software  process  will  run  on  the  CPU  to  copy  all  the  data  files  from  the 
disk  drives  to  the  tape  drive.  Once  the  transfer  to  tape  has  been  completed,  the  files  on 
the  disk  can  be  deleted,  freeing  up  disk  space  so  that  the  data  acquisition  process  can 
start  again.  Details  of  how  this  type  of  data  acquisition  system  could  be  applied  to  the 
AN/SQS-56  sonar  is  given  in  Appendix  D. 

Currently,  the  largest  storage  capacity  available  with  SCSI1  disk  drives  is  2.9  Gbytes. 
The  current  generation  of  workstations  typically  support  up  to  4  disk  drives,  giving 
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Figure  6:  Functional  block  diagram  of  a  computer-based  data  acquisition  system 

a  total  possible  storage  capacity  of  11.6  Gbytes  (if  such  a  computer  was  used  as  the 
host  CPU).  This  would  provide  sufficient  storage  capacity  for  more  than  3  hours  of 
data  at  a  data  rate  of  1  Mbyte/second,  and  more  than  12  hours  at  a  data  rate  of  250 
Kbytes/second.  However,  achieving  a  sustained  transfer  rate  of  1  Mbyte/second  in  this 
type  of  situation  is  unlikely  due  to  the  other  tasks  being  performed  by  the  CPU.  A  more 
realistically  achievable  data  rate  would  be  around  400-  500  Kbytes/second  sustained. 

At  present,  computer  tape  drives  have  not  reached  the  point  where  their  sustained 
transfer  rate  is  high  enough  to  make  them  usable  as  the  primary  storage  device  for  a 
computer-based  acquisition  system.  As  technology  improves,  this  situation  may  change, 
enabling  a  system  to  use  multiple  tape  drives  as  the  primary  storage  device.  Continuous 
storage  capability  would  then  be  achieved  by  changing  tapes  in  one  drive  while  the  other 
is  in  use. 

An  acquisition  system  of  this  type  would  be  very  flexible,  as  the  whole  process  is  under 
software  control.  It  would  simply  be  a  matter  of  re-configuring  the  software  to  cope 
with  different  sampling  rates,  number  of  channels  etc.  This  type  of  system  can  only 
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be  used  where  the  overall  data  rate  is  less  than  the  sustained  throughput  achievable 
with  current  computer  disk  drives.  As  such,  it  provides  an  effective  solution  to  data 
acquisition  problems  where  there  is  a  large  number  of  bandpass  signals  that  can  be 
sampled  at  a  rate  of  twice  their  bandwidth.  A  data  acquisition  system  similar  to  this  is 
being  developed  for  the  AN/SQS  56  sonar,  and  is  described  in  Appendix  D 


5  Summary 


An  overview  of  bandpass  signals,  and  special  techniques  for  sampling  bandpass  signals 
has  been  presented.  These  sampling  techniques  provide  an  efficient  method  of  sampling 
bandpass  signals  without  loss  of  information.  In  particular,  a  method  that  applies  to  a 
special  class  of  bandpass  signals  has  been  described  that  enables  a  bandpass  signal  to 
be  sampled  directly  at  a  rate  of  twice  its  bandwidth.  The  quadrature  representation  of 
the  bandpass  signal  can  be  extracted  from  the  sampled  data,  giving  a  different  method 
of  complex  demodulation  than  that  normally  used. 

By  using  this  method  of  complex  demodulation,  the  data  can  be  stored  immediately  after 
sampling.  The  lower  data  rate  achieved  using  this  sampling  method  enables  computer 
disk  drives  to  be  used  for  storage  in  place  of  the  large  magnetic  tapes  normally  used. 
Being  computer  based,  the  data  acquisition  system  uses  computer  storage  media,  which 
offer  a  number  of  advantages  over  a  tape  based  system.  Magnetic  tape  based  systems 
typically  use  special  encoders  when  recording,  requiring  the  tapes  to  be  replayed  through 
decoders  to  obtain  the  data  in  a  usable  form.  In  contrast,  the  method  described  in  this 
report  is  more  efficient,  providing  the  data  in  an  immediately  usable  form.  The  result 
is  savings  in  equipment,  work  effort,  time,  and  hence  cost. 

Another  advantage  of  this  type  of  system  is  its  flexibility.  It  is  easily  changed  to  suit 
different  data  acquisition  problems  (fewer  channels  at  a  higher  date  rate,  or  more  at  a 
lower  rate)  by  re-configuring  the  software.  The  system  is  also  easily  upgradable,  enabling 
the  capacity  of  the  system  to  be  increased  as  faster  and  larger  capacity  disk  drives  become 
available.  A  data  acquisition  system  of  this  type,  making  use  of  the  concepts  presented 
in  this  report,  is  being  developed  for  the  AN/SQS-56  sonar. 
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A  Complex  Demodulation 

A.l  Spectrum  of  a  complex  demodulated  signal 


Let  x(<)  be  a  bandpass  signal,  with  centre  frequency  /0,  bandwidth  \V,  and  spectrum 

*(/)• 


x(t)  =  xj(t)  cos(2x  fat)  -  x<j(f)  sin(27r/£,<) 
where  x/(f)  and  XQ(t)  are  obtained  as  follows 


( A- 1 ) 


x/(t)  =  2cos(2  rf0t)x(t) 
XQ(t)  =  2sin(2  irf0t)x(t) 


ip 


IP 


(A-2) 

(A-3) 


where  \lp  denotes  a  lowpass  filtering  operation. 

Define  F+(f  —  f0)  and  F_(/  +  /Q)  to  be  the  positive  and  negative  halves  of  the  spectrum 
of  the  bandpass  signal  x(f)  respectively.  Let  «=>  denote  the  Fourier  transform  operation, 
then 


i 


♦) 


« 


* 
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x(<)  X{f) 

X(f)  =  F4f-fo)  +  F.(f  +  fo)  (A-4) 

Also 

2cos(2t rf0t)  <*  S(f  -  fe)  +  6(f  +  /„)  (A-5) 

2sin(2x/0<)  «■  j£(/  -  /,)  +  j6(f  +  f0)  ( A-6) 

The  spectrum  of  the  in-phase  signal  x/(t)  is  given  by  the  convolution  of  (A-4)  and  (A-5), 
followed  by  a  lowpass  filter,  as  shown  by  (A-7)  below 


i 


LP 


*,(/)  =  *(/)*[«(/ -/.)  +  *(/  +  /.)] 

=  [F+(f  -  f„)  +  F-(f  +  fo)}  *  [S(f  -  fo)  +  6(f  +  fo)} 

=  F+(f  -  2 f0)  +  F+(f)  +  F-(f)  +  F-(f  +  2 fo) 

=  F+U)  +  F-(f) 


LP 


LP 


(A-7) 
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The  spectrum  of  the  quadrature  signal  xq(t)  is  given  by  the  convolution  of  (A-4)  and 
(A-6),  followed  by  a  lowpass  filter.  The  result  is  given  by  (A-8)  below. 


XqU)  =  XU)  *  hHf  -  fo)  +  jS(f  +  f0)i 


I LP 


=  [F+U  -  fo)  +  F-U  +  fo)}  *  hHf  -  fo)  +  j6(f  +  fo)} 


LP 


=  -jF+U  -  2 fo)  +  jF+U)  -  JF-U)  +  JF-U  +  -  fo) 

—  jF+{f)  ~  jF-{f) 


(A-8) 
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B  Overview  of  Quadrature  Beamforming 

B.l  Review  of  the  Conventional  Digital  Beamformer  Struc¬ 
ture 


Consider  the  conventional  digital  beamformer  shown  in  figure  B-l. 


y(m&) 


Figure  B-l:  Conventional  Time  Domain  Delay-Sum  Digital  Beamformer 

In  this  type  of  structure,  the  sampling  interval  6  (=  is  sufficiently  small  to  achieve 
the  delay  quantisation  required  for  beam  steering.  The  sampling  interval  to  satisfy  the 
Nyquist  criterion  is  A(=  which  is  often  an  integer  multiple  of  6  (this  is  not  a  general 
requirement  however),  i.e. 

A  ~  LS,  where  L  =  1,2, 3,4, .. .  (B-l) 

Referring  to  figure  B-l  for  simplicity,  assume  that  the  beam  output  y  is  produced  at  a 
sample  rate  of  A,  which  is  sufficient  for  waveform  reconstruction.  Note  that  the  beam 
output  can  be  computed  at  a  rate  K6,  where  K  is  an  integer  not  equal  to  L.  For  the 
former  case,  the  beam  output  y(mA)  is  found  from  the  weighted  sum  of  delayed  versions 
of  the  k  input  sequences.  The  delays  for  each  sensor  input  are  an  integer  multiple  of  the 
sampling  interval  6,  i.e. 

r,  =  M,6  (B-2) 

where  M;  is  an  integer  such  that  the  product  Mt6  gives  the  closest  approximation  to  the 
exact  delay  required  for  the  desired  beam  steering  angle.  For  the  case  where  the  beam 
output  rate  is  an  integer  sub-multiple  of  the  input  data  rate,  the  beamformed  output  is 
given  by  equation  (B-3) 
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(B-3) 


•  : 


k 

y(mA)  =  ^  a,Xi(mA  -  r.) 

i=i 

1=1 


where  a,  denotes  the  weighting  or  shading  applied  to  each  input  sequence.  The  total 
input  data  rate  for  this  structure  is  kf$,  and  the  output  data  rate  per  beam  is 
Implementation  of  equation  B-3  requires  real  additions  at  a  rate  of  (k  —  1  )/A  and  real 
multiplications  at  a  rate  of  kf&  per  beam. 

This  type  of  beamforming  structure  requires  high  input  data  rates,  increasing  the  com¬ 
plexity  of  the  hardware  required.  Interpolation  beamformers  [2]  result  in  reduced  hard¬ 
ware  complexity,  but  increase  the  computational  complexity.  However  for  bandpass 
signals,  even  interpolation  beamformers  use  sampling  rates  that  are  consistent  with  the 
highest  frequency  of  interest,  rather  than  consistent  with  the  bandwidth  of  bandpass 
signal.  The  following  section  describes  a  different  beamformer  structure  that  cam  be 
used  for  bandpass  signals. 
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B.2  The  Quadrature  Beamformer  Structure 


The  representation  of  a  bandpass  signal  in  terms  of  sampled  versions  of  the  in-phase  and 
quadrature  components  is  an  efficient  means  (in  terms  of  data  rate)  of  characterising 
the  signal.  A  quadrature  beamformer  is  a  beam  forming  structure  that  operates  on  this 
efficient  representation  of  a  bandpass  signal,  producing  beam  outputs  that  are  discrete- 
time  versions  of  the  in-phase  and  quadrature  components  of  the  beamformed  signal. 

If  the  beamformer  inputs  x i(n6)  are  bandpass  centred  on  frequency  /<>,  the  beamformed 
output  y(mA)  will  also  be  bandpass  centred  on  frequency  fa.  y{mA)  cm  be  expressed 
in  terms  of  its  in-phase  and  quadrature  components,  yi(mA)  and  yQ(mA),  as  shown  by 
equation  (B-4) 


•  •  4 


y(mA)  ~  yi(mA)  cos(2 ir  f0mA)  —  sin(2tr  f„mA)  (B-4) 

If  equation  (B-3)  is  re-written  so  that  the  sensor  inputs  are  expressed  in  terms  of  their 
in-phase  and  quadrature  components,  i.e. 
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1=1 


*  r 

y(mA)  =  a,  i/,(mA  -  r.)  cos(‘25r/0(mA  -  r,)J 

—  xg,(mA  -  r,)sin(2x/0(mA  -  r,)) 


»5> 


1=1 


i;,(mA  -  r,)  cos(27r/0mA)  cos(2ff/0r,) 

4-  ij,(mA  -  t,)  sin(27r /,mA)  sin(2r/0r, ) 

-  x^fmA  -  r,)sin(27r/0mAj)cos(25r/0rl) 

•f  xg,(mA  —  r.)  cos(2x/0mA)  sin(2r /0r, 


then  equating  equations  (B-4)  and  (B-5)  give 


yi(mA)  =  jx/,(mA  “  *») cos(2x/0r,)  +  xg,(mA  -  r,)sin(2r/0r,j 


(B-5) 


(B-6) 


and 


yg(mA)  =  Ea,jxg.(mA  -  r.) cos(27r/0r,)  -  x/,(mA  -  T,)sm{2xf0T,) 


(B-7) 


Equations  (B-6)  amd  (B-7)  show  that  the  in-phase  and  quadrature  components  of  the 
beam  output  can  be  obtained  from  the  in-phase  and  quadrature  components  of  the 
sensor  data. 

It  is  important  to  note  that  for  an  original  spectral  bandwidth  of  W,  both  yj  and  yg 
are  base  band  signals  of  bandwidth  -y.  Consequently,  the  sample  spacing  required  to 

adequately  represent  yj  and  yg  is  A  =  W~x,  rather  than  A  =  (2 (f0  +  y))  ,  (here  f0 

is  the  centre  frequency  of  the  bandpass  signal  as  before),  which  is  the  case  for  the  beam 
output  y. 

Hence  for  equations  (B-6)  and  (B-7),  the  sample  spacing  A  is  given  by  equation  (B-8) 


A  = 


W 


(B-8) 


As  before,  the  delays  r,  required  for  beam  steering  are  approximated  by  M,6,  where 
6  <  A.  However  because  the  in-phase  and  quadrature  components  of  the  sensor  data 
are  less  rapidly  varying  functions  of  time  them  the  sensor  data  itself  (lowpass  instead 
of  bandpass),  they  can  be  aligned  more  coarsely  in  time  prior  to  beamforming.  The 
sampling  interval  6  required  for  beamforming  needs  to  be  such  that  the  quadrature 
input  data  can  be  aligned  to  a  fraction  of  its  highest  wavelength  (y ).  Thus  6  depends 
on  W  rather  than  f0  +  y.  Reference  [3]  gives  a  comparison  of  this  type  of  approach 
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with  a  conventional  beamformer.  showing  the  effect  of  phase  errors  introduced  by  the 
coarser  delay  quantisation. 

The  information  regarding  the  carrier  frequency  f0  is  contained  in  the  sine/cosine  shading 
coefficients  that  weight  the  in-phase  and  quadrature  components  of  the  sensor  data  prior 
to  beamforming.  The  delay  r,  in  these  coefficients  must  be  specified  accurately,  but  this 
does  not  affect  A  or  6. 

As  the  in-phase  and  quadrature  components  of  the  sensor  data  are  both  lowpass  of 
bandwidth  they  can  be  adequately  represented  with  a  sample  spacing  of  A  given  by 
equation  (B-8),  with  the  vernier  delays  6  required  for  beam  forming  obtained  by  interpo¬ 
lation  [1],  [2].  This  results  in  both  the  input  and  output  data  rates  being  reduced  to  the 
minimum  allowed,  but  at  the  expense  of  increased  computation.  For  an  interpolation 
beamformer,  the  interpolation  can  be  carried  out  either  before  or  after  the  beamforming 
operation.  Figure  B-2  shows  a  quadrature  beamformer  with  the  interpolation  performed 
prior  to  the  beamforming. 


Figure  B-2:  Quadrature  beamformer  with  Pre-beamforming  interpolation 
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C  The  AN/SQS-56  Sonar  System 

C.l  Overview 

The  AN/SQS-56  sonar  system  is  a  member  of  the  DEI  160  series  of  surface  ship  sonars 
developed  by  Raytheon.  It  is  a  hull  mounted  sonar  that  can  operate  in  both  active  and 
passive  modes,  with  capabilities  for  detection,  classification,  localisation  and  tracking. 
When  operating  in  the  active  mode,  it  forms  36  equi-spaced  beams  covering  the  full  360° 
arc  around  the  ship. 

The  sonar  set  comprises  the  following  subsystems: 


•  Transducer  Array  Subsystem  which  provides  the  conversion  between  electrical  and 
acoustic  energy  and  the  physical  interface  between  the  system  and  the  medium. 

•  Array  Interface  Subsystem  which  provides  aperture  orders  to  the  transducer  and 
pre-amplifies  the  active  and  passive  signals  from  the  transducer. 

•  Transmitter  Subsystem  which  provides  the  transmitted  power  in  the  required  phase 
and  amplitude. 

•  Receiver  Subsystem  which  forms  the  received  active  and  passive  inputs  into  the 
required  beams,  suitably  processes  these  signals  and  delivers  digital  active  and 
passive  video,  as  well  as  audio,  to  the  console. 

•  Controller  Subsystem  whose  integral  general  purpose  computer  controls  the  func¬ 
tioning  of  the  system  in  response  to  stimuli  provided  by  the  operator  console.  It 
also  distributes  system  power,  and  provides  interface  to  certain  other  ship  systems. 

•  Control  Indicator  Subsystem  which  displays  the  processed  received  signals  and 
the  formatted  status  and  operation  data,  and  provides  the  interface  between  the 
operator  and  the  system  (i.e  the  console). 

•  Fault  Detection  and  Fault  Location  Subsystem  which  automatically  detects  and 
reports,  on  line,  system  faults  to  the  operator  at  the  console,  automatically  provides 
indications  of  possible  faults  by  audible  and  visible  alarms  and  localizes  these  to 
limited  areas  of  the  system. 

A  functional  block  diagram  showing  these  subsystems  (except  for  the  Fault  Detection 
and  Fault  Location  subsystem,  which  is  integrated  within  the  other  subsystems)  and 
their  major  inter-connections  is  shown  in  Figure  C-l. 

As  the  intended  application  of  the  type  of  data  acquisition  system  described  in  this 
report  is  to  the  active  data  stream  of  the  AN /SQS-56  sonar,  only  details  relevant  to  the 
active  mode  of  operation  will  be  considered.  When  operating  in  the  active  mode,  the 
sonar  transmits  on  one  of  three  nominal  centre  frequencies,  6.7,  7.5  and  8.4  kHz,  with  a 
bandwidth  of  600  Hz.  Further  details  of  the  operation  of  the  sonar  set  can  be  found  in 
(1)  -  [2]. 
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Figure  C-l:  AN/SQS-56  Sonar  -  Functional  block  diagram 


Figure  C-2:  Functional  block  diagram  of  Array  Interface  and  Receiver  Subsystems 
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Consider  the  functioned  block  diagram  of  the  Array  Interface  and  Receiver  Subsystems 
shown  in  figure  C-2. 

This  diagram  shows  that  the  active  receive  signals  from  the  hydrophones  which  make 
up  the  transducer  array  are  first  buffered,  then  summed  to  give  a  single  signal  per 
stave.  The  stave  sums  are  then  passed  through  a  bandpass  filter  (whose  centre  fre¬ 
quency  matches  the  selected  transmit  frequency),  followed  by  a  pre-amplifier/buffer, 
before  passing  through  isolation  transformers.  The  balanced  outputs  from  the  isolation 
transformers  pass  from  the  array  interface  subsystem  to  the  receiver  subsystems  via  a 
connecting  cable,  where  they  are  converted  to  digital  signals  using  a  geometric  analog  to 
digital  converter.  These  digital  signals  are  then  beamformed  via  a  digital  beamformer 
before  undergoing  further  processing. 

The  output  of  the  beamformer  would  seem  the  logical  place  from  which  to  extract  and 
record  the  data.  However,  to  allow  a  wider  range  of  processing  options  to  be  explored 
when  the  data  is  processed,  a  better  place  from  which  to  extract  and  record  the  data 
is  prior  to  beamforming.  Given  this,  the  best  place  for  the  data  to  be  extracted  is  at 
the  connection  between  the  array  interface  and  receiver  subsystems  (i.e.  after  the  active 
pre-amplifier  and  before  the  geometric  A/D  converter). 

Table  C-l  gives  the  characteristics  of  the  signal  conditioning  electronics  through  which 
the  active  data  signals  pass  within  the  array  interface  unit.  Further  information  can  be 
found  in  [1]  -  [2]. 


Characteristic 

Vedue 

Unit 

Remarks 

Voltage  Gain 

45.8/37.1 

dB 

Coherent /incoherent  gain/8  equal 
elements  referenced  to  element  input 

Beindwidth 

600 

Hz 

Meeisured  at  3  dB  down  points,  center 
frequencies  of  f\  =  6.7,  /2  =  7.5. 
fz  —  8.4  kHz.  Filters  shall  not  be  less 
than  50  dB  down  at  fx  ±  600  Hz. 

Dynamic  Range 

90 

dB 

Input  base  reference  level  =  -135dBV, 
1-Hz  band 

Output  Voltage 

+15 

dBV 

<  2%  distortion  into  1  K  ohm  paralleled 
with  1200  pF,  6-9  kHz 

Table  C-l:  Active  Pre- amplifier  Characteristics 

From  this  table  and  [1],  it  can  be  seen  that  the  36  active  signals  which  connect  the  array 
interface  to  the  active  receiver  have  the  following  characteristics 


•  bandpass,  centred  at  ftx  of  either  6, 7,  7.5,  or  8, 4  kHz,  with  3  dB  points  at  /tI±300 
Hz. 

•  signals  are  >  50  dB  down  at  ftx  ±  600  Hz 

•  transformer  isolated,  providing  balanced  signals 

•  maximum  signed  level  is  at  +15  dB  re  1  V  rms 


♦ 
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As  these  signals  at  the  desired  acquisition  point  are  bandpass  in  nature,  it  is  possible 
to  use  a  bandpass  sampling  technique  rather  than  conventional  sampling  techniques, 
resulting  in  a  lower  overall  data  rate  and  a  reduction  in  storage  requirements. 


♦ 
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D  Data  Acquisition  for  the  AN/SQS-56  Sonar 


D.l  Proposed  Sampling  Method  for  Active  Data  Stream 

In  Appendix  C,  the  point  identified  as  being  most  suitable  for  acquiring  the  data  was 
between  the  array  interface  and  receiver  subsystems.  However  at  this  point,  the  36 
active  signals  are  the  transducer  stave  sums.  Hence  the  sampling  method  used  must  be 
such  that  these  signals  can  be  passed  through  a  digital  beamformer  to  form  36  beams 
corresponding  to  the  36  beams  formed  internally  within  the  receiver. 

As  the  signals  at  the  identified  “pick  off”  point  are  bandpass  in  nature,  the  technique  of 
integer  band  shifting  on  sampling  (as  detailed  in  section  3.1)  can  be  used  to  sample  the 
signals,  providing  the  conditions  outlined  in  section  3.1  are  satisfied.  If  this  is  possible, 
the  method  of  complex  demodulation  detailed  in  section  3.2  can  be  used  to  convert  these 
sampled  signals  to  sampled  representations  of  the  in-phase  and  quadrature  components 
of  the  “stave  sum”  sign jlIs.  These  quadrature  signals  can  then  be  processed  by  a  quadra¬ 
ture  beamformer  to  form  the  36  beams  required  to  emulate  the  36  active  beams  used  by 
the  sonar  set.  An  outline  of  quadrature  beamforming  is  given  in  Appendix  B. 

Consider  first  the  case  where  the  sonar  is  operating  on  its  highest  frequency,  ftz  =  8.4 
kHz.  The  stave  sum  signals  are  at  least  50  dB  down  at  8.4  ±  0.6  kHz.  Assuming 
a  minimum  signal  to  noise  ratio  of  50  dB  for  further  processing,  the  minimum  signal 
bandwidth  that  can  be  tolerated  for  sampling  will  be  1200  Hz  (/*,  =  7.8  and  fu  =  9.0 
kHz),  requiring  a  minimum  sample  rate  of  at  least  twice  this.  As  the  upper  and  lower 
frequency  bounds  fu  and  fi  actually  denote  the  stop  band  edges,  equation  (14)  of 
section  3.1  applies.  If  W  denotes  the  signal  bandwidth,  fi,  fu  and  f,  can  be  written  as 
follows 


Jl 

fu 

f. 


w 

-ft:  j 

(D-l) 

w 

-  /«x  +  y 

(D-2) 

=  2W 

(D-3) 

If  equation  (D-3)  is  re-written  to  express  W  in  terms  of  /,,  and  then  substituted  into 
equation  (D-2),  which  is  then  substituted  into  equation  (12)  of  section  3.1,  an  expression 
for  f3  in  terms  of  k  and  ftx  is  obtained. 


/,  =  JtX'i  where  k  —  1,2, 3,4, .. .  (D-4) 

k  ~  4 

Table  D-l  shows  /,,  /i,  fu  and  W  for  various  values  of  k 

It  was  mentioned  above  that  the  minimum  bandwidth  which  provides  adequate  attenua¬ 
tion  for  out  of  band  signals  is  1200  Hz,  and  so  only  the  first  3  solutions  given  in  table  D-l 


H 

Sampling 
Frequency,  f,  (Hz) 

Lower  Frequency 
Bound,  }i  (Hz) 

Upper  Frequency 
Bound,  fv  (Hz) 

Bandwidth 
W  (Hz) 

B 

11200 

5600 

11200 

MM 

B 

7200 

9600 

WBm 

2 

3055 

7636 

9164 

H 

2240 

7840 

8960 

In 

Table  D-l:  Possible  sampling  frequencies  for  integer  band  shifting,  fix  =  8400  Hz 

apply,  with  the  third  solution  providing  the  best  compromise  in  terms  of  data  rate  and 
out  of  band  attenuation. 

Similar  tables  for  the  two  other  transmission  frequencies  of  6.7  and  7.5  kHz  are  shown 
in  tables  D-2  and  D-3.  The  figures  given  in  these  tables  were  evaluated  using  similar 
method  to  that  used  for  table  D-l 


H 

Sampling 

Frequency,  /,  (Hz) 

Lower  Frequency 
Bound,  Jl  (Hz) 

Upper  Frequency 
Bound,  fv  (Hz) 

Bandwidth 
W  (Hz) 

i 

•si 

4467 

8933 

4467 

2 

5743 

7657 

1914 

3 

2436 

7309 

1218 

1787 

7147 

893 

Table  D-2:  Possible  sampling  frequencies  for  integer  band  shifting,  ftx  —  6700  Hz 


Sampling 
Frequency,  f,  (Hz) 

Lower  Frequency 
Bound,  ft  (Hz) 

Upper  Frequency 
Bound,  fv  (Hz) 

Bandwndth 
W  (Hz) 

1 

10000 

5000 

10000 

5000 

2 

4285 

6429 

8571 

2143 

3 

2727 

6818 

8182 

1364 

m 

2000 

7000 

8000 

1000 

Table  D-3:  Possible  sampling  frequencies  for  integer  band  shifting,  ftx  =  7500  Hz 

For  the  cases  outlined  in  tables  D-2  and  D-3,  only  the  first  3  solutions  satisfy  the  min¬ 
imum  bandwidth  requirement,  with  the  third  solution  giving  the  lowest  sampling  fre¬ 
quency. 

Note  that  tables  D-l  to  D-3  are  for  the  case  where  /,  =  kfv-  A  similar  set  of  tables  for 
the  other  case  of  f,  =  kfi  can  be  obtained  by  expressing  /,  as  a  function  of  f,z  and  k, 
and  working  through  the  numbers  (an  expression  is  obtained  similar  to  the  method  for 
equation  (D-4),  except  equation  (9)  from  section  3.1  is  used  in  place  of  equation  (12)  in 
the  substitution).  This  option  would  also  require  the  method  outlined  in  section  3.2  to 
be  modified  to  obtain  the  in-phase  and  quadrature  signals  as  outputs. 
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These  tables  show  that  for  each  of  the  3  transmit  frequencies,  integer  band  shifting  on 
sampling  is  possible  without  any  unintentional  aliasing.  This  means  that  the  method  of 
complex  demodulation  described  in  section  3.1  is  applicable  to  the  active  data  signals, 
resulting  in  a  much  lower  overall  data  rate  than  would  be  possible  using  conventional 
sampling  techniques.  For  all  three  transmit  frequencies,  the  higher  sampling  frequency 
will  give  the  best  signal  to  noise  ratio  (maximum  out-of-band  attenuation,  thus  mini¬ 
mizing  objectionable  aliasing),  but  result  in  higher  data  rates.  The  final  sampling  rate 
used  will  be  the  highest  rate  the  data  acquisition  system  can  cope  with. 


D.2  Data  Acquisition  System  Proposal 


The  overall  function  of  the  data  acquisition  system  is  to  digitize  the  active  data  from 
the  sonar  and  store  it  for  subsequent  analysis.  Because  of  the  nature  of  the  bandpass 
active  signals  on  the  AN/SQS-56  sonar,  it  is  possible  to  sample  the  signals  at  a  rate 
of  twice  their  bandwidth,  as  described  in  Appendix  D.l.  This  sampling  method  results 
in  data  rates  of  the  order  of  several  hundred  Kbytes  per  second  for  all  36  active  data 
signals.  Combining  this  sampling  method  with  the  type  of  system  architecture  described 
in  section  4  will  give  a  data  acquisition  system  whose  stored  output  is  computer-readable, 
and  ready  for  immediate  analysis. 

The  functional  requirements  of  a  data  acquisition  system  for  the  AN/SQS-56  sonar 
capable  of  the  above  are  briefly  as  follows 


•  minimum  of  12  bits  per  sample  (preferably  16  bits) 

•  synchronously  acquire  data  for  all  36  data  channels  at  a  rate  sufficient  for  beam¬ 
forming 

•  store  sever  ad  hours  of  acquired  data 

•  data  to  be  stored  in  computer- readable  form,  ready  for  immediate  analysis 

•  either  store  the  data  as  one  data  file  per  ping  or  mauntain  a  time  log  of  the  pings 
so  that  the  data  can  be  separated  on  a  per  ping  basis  at  some  later  stage. 


•  • 


A  system  consisting  of  the  following  components  would  be  able  to  satisfy  these  require¬ 
ments. 


•  5  ICS-140-8A  8  channel  VMEbus-  based  data  acquisition  cards 

•  1  SPARC  IE2  VMEbus  CPU  card  running  VxWorks3  operating  system  software  ( 

•  1  SPARCstation4 

2SPARC  is  a  registered  trademark  of  SPARC  International,  Inc. 

3VxWorks  is  a  trademark  of  Wind  River  Systems,  Inc. 

4SPARCstation  is  a  registered  trademark  of  SPARC  International,  Inc.,  licensed  exclusively  to  Sun 
Microsystems,  Inc.  •  < 
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•  1  or  more  large  capacity  SCSI  disk  drive  (>  1.2  Gbytes) 

•  8mm  (Exabyte4)  tape  drive 

A  functional  diagram  showing  how  these  components  inter-connect  to  form  the  data 
acquisition  system  is  given  in  figure  D- 1 . 


Active  Data  Active  Data  Active  Data  Active  Data  Active  Dma 

Staves  1  -  8  Staves  9  - 16  Staves  17  -  24  Staves  25  -  32  Staves  33  -  36 


Ping  Signal 


Figure  D-l:  Proposed  SQS-56  Data  Acquisition  System 

The  5  ICS-140-8A  data  acquisition  cards  and  the  SPARC  IE  CPU  card  are  housed  in  a 
single  VME  chassis.  An  ethemet  cable  connects  the  SPARC  IE  CPU  to  the  SPARCsta- 
tion,  which  has  a  large  capacity  disk  drive  and  8mm  tape  drive  as  storage  peripherals. 
The  SPARC  IE  will  be  running  VxWorks  software,  a  reed  time  UNIX6-compatible  op¬ 
erating  system,  while  the  SPARCstation  will  be  running  standard  SunOS7. 

5  Exabyte  is  a  registered  trademark  of  EXABYTE  Corporation 
®UNIX  is  a  registered  trademark  of  AT&T  Bell  Laboratories 
7SunOS  is  a  trademark  of  Sun  Microsystems  Inc 


The  36  stave  sum  signals  are  extracted  from  the  interconnection  between  the  array  inter¬ 
face  and  receiver  subsystems,  and  then  passed  through  some  simple  signal  conditioning 
electronics.  This  signal  conditioning  is  mainly  to  attenuate  the  stave  sum  signals  so  that 
the  signal  level  is  compatible  with  the  analog  inputs  on  the  ICS-140-8A  data  acquisition 
cards.  Following  this,  the  36  ‘"conditioned  signals  are  connected  to  36  of  the  40  analog 
inputs  on  the  ICS- 140-8 A  cards  (all  8  inputs  on  the  first  4  cards  will  be  used,  and  only 
4  on  the  fifth  card). 

The  SPARC  IE  CPU  will  be  running  Vx  Works,  under  which  will  be  running  a  number 
of  software  processes  to  control  the  operation  of  the  ICS  cards  and  the  data  acquisition 
process  itself.  Prior  to  any  acquisition  taking  place,  each  card  will  be  configured  ready 
for  acquisition.  Once  configured,  the  acquisition  process  can  be  started  at  any  time.  The 
process  that  controls  the  data  acquisition,  reads  the  data  acquired  by  the  ICS  boards 
over  the  VME  bus,  and  then  writes  this  data  to  a  disk  file  on  the  SPARCstation  via  the 
ethernet  connection.  When  a  start  of  ping  signal  is  received  by  this  process,  the  current 
data  file  is  closed  and  a  new  one  opened.  Thus  there  will  be  a  unique  data  file  for  each 
ping,  containing  all  the  data  samples  from  all  36  transducer  sta\es  :  r  that  ping. 

At  the  end  of  the  trial  period,  all  the  data  files  stored  on  the  hard  disk  on  the  SPARC¬ 
station  (several  hundred  to  approximately  a  thousand  files),  will  be  copied  from  the 
disk  to  the  8  mm  tape  drive  for  storage.  Once  the  copying  to  tape  has  been  completed, 
the  files  can  be  deleted  to  free  up  disk  space  in  preparation  for  the  next  trial.  As  the 
data  files  stored  on  the  tapes  were  originally  written  by  a  Vx Works  process  they  will  be 
UNIX  files  and  thus  will  be  computer-readable  requiring  no  further  processing  before 
being  analysed. 

As  the  whole  acquisition  and  storage  process  will  be  under  software  control,  it  will  be 
a  simple  task  to  select  the  appropriate  sampling  frequency  (given  in  tables  D-l  -  D-3) 
depending  on  the  current  transmit  frequency.  This  flexibility  together  with  the  data 
storage  method  used  makes  a  system  such  as  that  outlined  above  an  ideal  solution  for 
acquiring  and  storing  the  active  data  stream  from  the  AN/SQS-56  sonar. 


D.3  Acquisition  and  Storage  of  Ancillary  Ship  Data  • 

To  enable  any  analysis  carried  out  on  the  recorded  trials  data  to  be  interpreted  correctly. 

knowledge  is  required  of  the  ships  behaviour  throughout  the  period  of  the  trial  (i.e. 

speed,  heading,  roll  etc).  This  information  is  available  on  the  ship  from  the  Naval 

Combat  Data  System  (NCDS).  A  system  that  can  decode  the  information  recorded  * 

by  the  NCDS  is  already  available.  This  system  extracts  information  from  the  NCDS 

recordings  and  stores  it  on  a  floppy  disk  ready  for  analysis  at  some  later  stage.  To  use 

this  system  all  that  would  be  required  is  to  synchronise  the  times  stored  by  the  NCDS 

with  those  times  stored  by  the  data  acquisition  equipment.  This  would  enable  the  ship’s 

behaviour  during  each  individual  ping  to  be  known  accurately.  * 
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The  function  of  a  data  acquisition  system  is  to  digitise  data  from  a  number  of  separate  channels,  and  store  these 

samples  for  later  analysis.  Sonar  systems  are  typically  multi-channel,  and  the  associated  data  acquisition  0  0 

systems  are  required  to  handle  data  at  high  data  rates.  The  high  data  rates  involved  in  digitising  bandpass 

sonar  data  usually  mean  that  the  only  viable  storage  methods  available  are  large  magnetic  tapes,  with  data 

recorded  on  high-speed  tape  recorders.  One  of  several  disadvantages  of  this  type  of  storage  method  is  that  the 

types  have  to  be  replayed  through  special  decoders  to  obtain  the  data  in  a  form  ready  for  analysis.  This  is  a 

time  consuming  process,  delaying  the  analysis  and  interpretation  of  die  data.  An  alternative  method  of 

sampling  bandpass  sensor  data  which  results  in  a  much  lower  data  rate  is  presented  in  this  report  The  •  • 

resulting  lower  data  rate  allows  advantage  to  be  taken  of  the  recent  improvements  in  computer  storage 

technology,  enabling  a  data  acquisition  system  to  be  constructed  that  is  computer-based,  and  utilises  computer 

storage  devices  as  its  primary  storage  media.  The  advantage  of  this  storage  method  is  that  the  data  is  in 

computer  readable  form,  ready  for  immediate  analysis.  This  work  was  carried  out  in  support  of  Navy  Task 

NAV  92/401  and  an  application  of  the  theory  of  the  AN /SQS  56  sonar  is  presented  in  the  Appendices.  _  _ 
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